Covariant nucleon wave function with S, D, and P-state components 
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Expressions for the nucleon wave functions in the covariant spectator theory (CST) are derived. 
The nucleon is described as a system with a off-mass-shell constituent quark, free to interact with 
an external probe, and two spectator constituent quarks on their mass shell. Integrating over the 
internal momentum of the on-mass-shell quark pair allows us to derive an effective nucleon wave 
function that can be written only in terms of the quark and diquark (quark-pair) variables. The 
derived nucleon wave function includes contributions from S, P and D- waves. 



I. INTRODUCTION 

In this work we use the covariant spectator theory 
(CST) P~H5] to determine the structure of the valence 
quark contributions to the wave functions of the nucleon. 
In the CST, baryon systems consist of an off mass-shell 
constituent quark, free to interact with electromagnetic 
fields or other probes, and two noninteracting on mass 
shell constituent quarks that are spectators to the in- 
teraction. Since the interaction does not depend on the 
internal momentum of these on-shcll spectators, we can 
integrate over their internal momentum and express the 
effective matrix element in terms of a nucleon composed 
of an off-shell quark and an on-shell quark pair (or di- 
quark) with an average mass m s , which becomes a pa- 
rameter in the wave function. 

Previously we have assumed a pure S-wave structure 
for the wave functions; here we add contributions asso- 
ciated with D and P-wave components. We begin by 
discussing the general form of the CST matrix elements. 
Then, for each angular momentum component, nonrel- 
ativistic wave functions are constructed first, and then 
generalized to relativistic form. 



II. CST MATRIX ELEMENTS 

A. Relativistic impulse approximation 

In the relativistric impulse approximation (RIA), it is 
assumed that the interaction is well described by a sin- 
gle quark operator O a , and that interactions involving 
a pair of quarks (i.e. exchange or interaction currents) 



can be neglected. In this case, the baryonic matrix ele- 
ments in the CST can be written (for a brief discussion 
of corrections to the RIA, see Ref. [B]) 

i=l \j\ e Jk i k e 

where the first sum is over all three possible choices for 
the interacting quark, and the three-momentum integra- 
tions and helicity sums are over the momenta and helic- 
ities of the on-shell spectator quarks, with i, j, I in cyclic 
order, so that, for example, if the third quark is interact- 
ing, then the (12) pair are the spectators. The matrix 
element is 

(O^^^x^x+iP+^jk^O^x^xJP-^jh) (2.2) 

with P_ (P+) and A_ (A+) the four-momenta and helic- 
ity of the incoming (outgoing) baryon, and ^ the baryon 
wave function. We have suppressed the Dirac indices of 
the off-shell quark. The matrix element is illustrated in 
Fig. [T] In the CST the spectator quarks are constrained 
to their positive energy mass-shell, so the covariant vol- 
ume integral is (with j '• = 1 and £ = 2) 

j ^ j d^k^,^ _ kl)6+[ml _ ^ 

J (2n)HE 1 E 2 [ ■ ' 

where the four-momenta of the on-shell quarks are 
h = {E 1 ,k 1 } 

k 2 = {E 2 ,k 2 }, (2.4) 
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FIG. 1. (Color on line) Diagrammatic representation of the RIA 
approximation to the CST matrix element of the operator O m . The 
quarks with four moments k\ and &2 are on-shell spectators (rep- 
resented by the X ) . 



with Ei = y'mf+kf'. 

In this paper we will first assume the dressed masses of 
the three quarks are equal, so that rri\ = m 2 = m 3 = m. 
(Later we will consider the case when m u 7^ m^-) Since 
the color factor (suppressed) is fully antisymmetric, when 
the particles are identical this means that the remaining 
matrix element must be symmetric under the interchange 
of the three quarks. In this case the three terms for 
different i are identical, and we may write the full result 
as three times the result for the (a rbitr ary) choice of k\ 
and fc 2 as spectators transforming (12. II) to 



Aj 



Ai A 2 Jklk: 



(2.5) 



This means that it is only necessary for the wave function 
*&(P, ki,k 2 ) to be symmetric under the interchange of 
quarks 1 and 2; symmetry under the interchanges of the 
other quarks, 1 <4 3 or 2 <4 3, was used when the full 



result (2.1) is simplified to (2.5). 



B. 



Relativistic definition of the quark-diquark 
wave function 



As shown in Fig. [T] the dependence of the matrix 
element on the relative momentum of the two on-shell 
quarks is determined by the wave function only, and not 
by the structure of the operator , which depends only 
on the momenta k%± of the initial and final off-shell quark 



k 



3± = 



(2.6) 



Since the wave function will be determined phenomeno- 
logically, very little is lost by averaging over the rela- 
tive momenta of the two on-shell quarks and introduc- 
ing a new wave function that depends only on the total 
four-momentum of the on-shell quark pair (which will be 
called a diquark). In order to do this we introduce the 
diquark momentum variables 



k 2 ) 



(2.7) 



which, together with the fixed total momentum P = k\- 



these would be the Jacobi coordinates.) The diquark 
four-momentum k = k\ + fc 2 has a mass s = (k\ + fc 2 ) 2 \ 
which is initially unconstrained. In this notation, the 
integrals ( |2.3| ) can be re-expressed as integrals over k, s, 
and the direction of the internal relative three momentum 
r = ^(ki — k 2 ) of the diquark 



s k 



4(2tt) 3 



ds 



4m; 



4m 2 



d 3 k 



(27r)32£ s ' 



(2. 



Is 



where m q is the dressed quark mass, and E s = v- 



s + k 2 

is the energy of a diquark of mass y/s, and the angular 
integrals d£l T are written in the rest frame of the diquark. 
The relation (2.8) is discussed in Appendix |A} 



Using the relation (2.8) the matrix element (2.5 1 can 
be written 



O 



A + A_ 



(P+,P-) 



A 



*AxJP+,k)0^ AX _(P_,k) (2.9) 



where the sum over A includes all possible polarization 
pairs Ai, A 2 of the diquark, and the new density in (2.9) 
is related to the density in (2.5) by 



Ai A 2 



*A 1 A 2 ,A+0P+, k X k 2 ) ® *A!A 2 ;A_ (P-, kxk 2 ) 



E 

A 



*AA + CP+,fc)®*AA_OP-,fc) 



,(2.10) 



where the operator Cgi shows where the Dirac operator O m 
is to be inserted. This equation gives the precise relation- 
ship between the quark-diquark wave function, denoted 
by 'I'a.aC-Pj ^)j an d t ne full three-quark wave function, 
^\i\ 2 ,x{P , k±k 2 )- (Note that these two wave functions 
are distinguished from each other only by their list of ar- 
guments.) Equation ( |2.10 1 shows that the quark-diquark 
wave function is obtained from the full three-quark wave 
function by averaging over the directions of the relative 
three momentum r (in the rest system) and replacing 
the integral over the continuous diquark mass, <Js, by 
its mean value, m s , which now becomes a parameter of 
the theory. All remaining factors from the integral J s 
are absorbed into the normalization of the wave func- 
tion. After the average has been carried out, the diquark 
energy becomes E s = \fm1 + k 2 . 



We emphasize that, in the RIA, the replacement (2.10) 
is exact, as long as the effective diquark mass m 



is 



treated as a parameter. Later we will constrain the values 
of this parameter, and in so doing we make an approxi- 
mation. 



Nonrelativistic definition of the quark-diquark 
wave function 



The nonrelativistic limit of (2.8) is obtained by assum 



&3, are a complete set. (If k — > k± + fc 2 — 4P, ing the quark and diquark masses are very large, and 
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absorbing the mass factors into the wave function nor- 
malization. In this case 



k 2 ) k 2 



4m„ 



Ay 2 



(2-11) 



For very large masses (E s —> m s ; s — > 4m^), we can write 



s k 



1 r d 3 k r dfif 

16 J (2tt)3 J (2tt)3 
d 3 k f d 3 r 



Am 2 



(2tt) 3 J (2tt) 



(2-12) 



where, in the last step, the integral over s has been re- 
placed by an integral over the magnitude of the relative 



momentum, r = r 



, using (2.111 



The result (2.12| also follows directly from the trans- 



formation from ki,k2 to k,r, which gives 



d 3 k x d 3 k 2 
(2tt) 6 



d 3 k 



(2tt)3 J (2n) 3 ~ J k 



NR 



(2.13) 



Hence the nonrelativistic definition of the diquark wave 
function is identical to (2.10), with 



d 3 r 
(27P 



(2.14) 



Since the details of how the averaging over the internal 
diquark structure depends on the individual angular mo- 
mentum components of the wave function, this discussion 
is postponed until Sec. |IV| 



III. ELECTROMAGNETIC INTERACTION 
AND NORMALIZATION OF THE WAVE 
FUNCTION 

The interaction of the photon with a constituent quark 
is decomposed into Dirac and Pauli components. 



r(Q 2 ) = 3i r 



■32- 



2M 



(3.1) 



where M is the nucleon mass. Note that the constituent 
quarks can have anomalous magnetic moments, and that 
we add the term —fq^/q 2 to insure that the current is 
always conserved [7j. In this paper we will limit discus- 
sion to u (I z = t = +|) and d (t = — i) quarks only, 
in which case the quark current j 7 ' and the quark form 
factors ji are operators in isospin space. The Dirac (ji) 
and Pauli (j'2) form factors are therefore decomposed into 
two independent isospin structures 



3 l {Q 2 ) = lh+{Q 2 ) + \f l -{Q 2 )T^ 



(3.2) 



where T3 is the ^-projection of the isospin operator. This 
equation defines the isoscalar and isovector form factors 
that can also be expressed as the u and d quark form 



factors. In the CST quark model the functions fi±(Q 2 ) 
are parametrized using a vector meson dominance repre- 
sentation [2], 



Using (2.9) and the operator (3.1), the baryon form 



factors become 



= 3 E / ^AX+iP+WiQ 2 )*^-,®, (3.3) 

A Jk 

For the nucleon, at Q 2 = 0, when P + = P_ = (M, 0, 0, 0) 
and 

/(0) = eo(i + ir 3 ) 7 = J1 (0) 7 ° (3.4) 

(because of the presence of the factor eo, this operator at 
Q 2 = gives the renormalized quark charge), we require 
that the wave function be diagonal in the nucleon polar- 
ization and normalized to the correct nucleon charge: 



J° x+x _(P 1 P)^Q = ±(1 + t 3 )5 x 



(3.5) 



It is always possible to normalize the wave functions so 
that the proton charge is unity, but to obtain zero for the 
neutron charge places constraints on the structure of the 
wave function. 



IV. WAVE FUNCTIONS OF THE NUCLEON 

A. Total relativistic wave function 

The total relativistic quark-diquark wave function will 
be written as a sum of S, P, and D-wave components 

= n s * s AX (P, k) + n P ^ x (P, k) + n D *% x (P, k) (4.1) 

where *Bj[ x (P, k) is the L=S, P, or D-state wave function 
to be defined in the next subsections. 

We will require that each of these components be sep- 



arately normalized to the nucleon charge, Eq. (3.5). This 



requirement will be decisive in fixing the structure of the 
wave functions. When each component has been sep- 
arately normalized, the S-state normalization constant, 
?t-5, will be determined by the P-state mixing parameter, 
np and the D-state mixing parameter no 



n s 



1 



(4.2) 



B. S-state component 

1. Spin-isospin part of the nonrelativistic wave function 

The construction of the nonrelativistic S-state proton 
wave function was given previously [5] , but here the con- 
struction will done differently. We obtain the same result 
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as we did before, but this way of doing things will be 
helpful when we extend the discussion to the D-state. 

The wave function is written as a product of a momen- 
tum space wave function symmetric in k 1 and k 2 (this is 
the only possibility for an S-state) and a spin-flavor part 
also symmetric under the interchange of quarks 1 and 2. 
There are two possible structures that contribute to the 
symmetric spin-flavor part: the (0,0) component which 
is a direct product of an operator antisymmetric in the 
flavor (isospin 0) and antisymmetric in the spin (spin 0), 
and the (1,1) component symmetric in the flavor (isospin 
1) and symmetric in the spin (spin 1). The isospin and 
1 components can be written in the following form 



i 



with x the two-component isospinor of the nucleon 



(4.3) 



+ A 
x +2 = 



X 2 = 



(4.4) 



and ^£ the isospin vector of the isospin-one diquark, de- 
fined in the usual way (in rectangular coordinates) 





(4.5) 



(This notation differs slightly from Ref. [2], but the re- 
sults are the same.) The spin one vectors are normalized 
to 



5> 



6e, 
8 a 



(4.6) 



where the ith component of the vector £| is denoted (&)». 
In the nonrelativistic limit, the spin wave functions and 
operators have precisely the same form. 

Using the operator notation (4.3) the total wave func- 



tion of the nucleon can now be written in a general form 
which displays its dependence on all three quark mo- 
menta, but still describes the flavor and spin of the three 
quarks in terms of the quark-diquark description already 
developed. By separating explicitly the spin isoscalar 

S 

diquark contribution, , from the spin 1 vector di- 

S 1 

quark contribution, , we have 



AM P; ki, k a ) = cos 6 S <j> u ¥£' U (P; k 1; k 2 ) 



-sin# s 1 *^(P;ki,k 2 ) 



(4.7) 



where A is the spin projection of the nucleon, A the spin 
projection of the spin-1 diquark, and ^/^(P; ki, ka) is 
the S-state spin-space part of the wave functions (with 
s the spin of the diquark and no A dependence in the 
spin-0 component). Explicitly 



^•"(Pjk^ka) 



ka) = -^s(P; ki, k 2 )a ■ e* A ||, A) , (4.8) 



where we have introduced our notation for the diquark 
spin functions, modeled after the isospin operators ( |4.3[ ), 
with e A the spin vector for the spin-1 diqu ark [with its 
components in spin-space defined as in (4.5 1], and ||,A) 
the initial nucleon spin state with polarization A. At this 
point we assume that the same momentum space wave 
function, $5, accompanies both of the spin parts of the 
wave function; later this assumption will be relaxed as 
the wave function is determined from phenomenological 
fits to the DIS data. We emphasize that, at this stage, 
the two components of the wave function have spins equal 
to their isospins, but it is best to label these components 
by their spin because later (in Sec. |V|) we will break the 
isospin invariance. 

Using this wave function, the nonrelativistic matrix 
element of the charge operator, (3.5) becomes 



Q=i(l + T 3 R +A _ 

={ cos 2 S j A + sin 2 6 S f}s x+>x _ e Q Af s (4.9) 

where the isospin matrix elements are 

^3(0^(0)0° = ^ + IN 

j 8 = 3]T (01)^1(0) 4>\ = J> • r) M0) (C e ■ r) 

i i 

= \e Q (1 - r 3 ) , (4.10) 



where ji(0) is the quark charge operator (3.4), the nor- 
malization integral is 



A/"c = 



d 3 k 



(2tt) 3 J (2tt) 



$1(0;^, k 2 ), (4.11) 



and we used the fact that the spin operators in (4.8) are 
rcnormalized. Choosing Ms = 1/eo, the charge operator 
(4.9) becomes 



Q= { cos 2 6 S 



1-3 



3 cos 2 





1 1 


+ sin 2 6 S 


2 ~ 2 T 3 


— sin 2 9s 





(4.12) 



Hence an equal mixture of diquark spin-0 and spin-1 com- 
ponents (6s = tt/4), which we have used in our previous 
work, is required by the demand that the neutron charge 
be zero. We will find this requirement useful in the con- 
struction of the D-state below. 



The integral (4.11 ) provides an explicit example of how 
the quark-diquark wave function emerges from an aver- 
age over the internal momenta of the diquark, as shown 
in general in Eq. ( 2.10 ). Since the momentum wave func- 



tion must be symmetric in k\ and fca, we may choose its 
argument to be 



Xr 



-(fc^rMki-iP^ + Cka-lP) 2 



= Ik 2 
= l(k- 



2r 



-fk-P 
) 2 +2r 2 



2p2 



(4.13) 



5 



With this choice, the quark-diquark wave function in an 
arbitrary frame is defined by the relation 



d 3 k 
(2tt) 3 



d 3 r 
(2tt) 3 
d 3 k 

d 3 k 
(2tt) 3 



®S [Xnr (k,P+,r)]$ 3 [Xnr (k,P-,r)] 
®S [Xnr (k,P+,r)]$ S [Xnr (k,P-,?)] 



MP+,k)0s(P-,k) 



(4.14) 



where, in the second step we replace r by and average 
value f, and the last step gives the precise relation be- 
tween the quark-diquark wave function, <j)s (P,k), and 
the three-quark wave function, $5(P;k 1 ,k 2 ), analogous 
to the definition given in Eq. ( 2.10 1. With this definition, 



the normalization of the quark-diquark wave function is 

/d 3 k 
„4(0,k) = l. (4.15) 

Note that the normalization of the wave function com- 
pensates for the renormalization of the quark charge, en- 
suring that the proton charge is correct. 



2. RelativisUc wave function 

The relativistic generalization of the nonrelativistic 
wave function is straightforward, and has been discussed 
extensively in Ref. [2J. The wave function we use is 



[<j>°u{P,X) — <p 1 (e%)*U a (P,Xj\ 



x MPk). 



(4.16) 



where the relativistic quark-diquark wave function, 
^^(P, k), is distinguished from the nonrelativistic, 
three-quark wave function, ^^(P; fc, r) only by its ar- 
guments, and U a is 



C/ Q (P,A) = - 7 = 75 



7a 



Pa 

M 



t (P,A), (4.17) 



with u(P, A) the free nucleon spinor with four-momentum 
P and spin projection A, £a is the polarization vector 
of the spin-1 diquark state with polarization A in the 
direction of P (and subject to the constraint s^P^ = 0), 
and the 0,1 are the flavor wave functions of the quarks 
in a (12)3 configuration with the (12) pair in an isospin 
zero or one state [see Eq. (4.3)]. The polarization vectors 
are discussed in detail in Refs. [S . The states have 



been normalized so the superposition (4.16) corresponds 
to an equal mixture of (0,0) and (1,1) components, as 
discussed above. 

The spatial part of the S-state wave function will be 
fixed by comparison with the deep inelastic scattering 
(DIS) data as discussed in the accompaning paper [3]. 
The model used in earlier work had the form 



iVo 



™ s (/3i+x)(/3 2 + x) 



(4.18) 



where 



_ (M - m s ) 2 - (P - k) 2 _ 2P ■ k 2 



Mm, 



Mm* 



(4.19) 



and /?2 > Pi > are range parameters with the normal- 
ization constant Nq chosen so that 



(4.20) 



with the covariant integration defined by Eq. (2.8) (with 
s = m 2 s ). Note that, in the nonrelativistic limit, x re ~ 
duces to 



M 2 



J-(k-*P 

Am 2 „ \ 3 



2m 2 



Mm. 



Xnr(k,P,0) 



k P 

m, M 



(4.21) 



where in the last step, to display the relation between 
X and Xnr: we made the approximation m s ~ 2m q and 
M ~ 3m q , true when all momenta (and binding energies) 
are much smaller that m q . For future reference we note 
now that 



(P-k)P 
P 2 



(P ■ fc) 2 
M 2 



(4.22) 



With the normalization (4.20) the nucleon charge is 



A 



Q = ■> > ; / vLiPVnioWfLiPk), (4.23) 

Ik 



where ji(0) is the quark charge operator defined in 
Eq. (|3"!4|). This agrees with Eq. (2.14) of Ref. [1^ (with 



the sum over the diquark polarization A extended to in- 
clude the spin 1 diquark with polarization A and the spin 
diquark as a separate term). 

We now discuss the construction of the D-state wave 
function. 



C. D-state component 

1. Three-quark nonrelativistic wave function 

In order to produce a spin 1/2 state through the cou- 
pling of a D-wave operator with a state of total quark 
spin S, the quark spin S must be 3/2. This is a purely 
symmetric spin state. Hence, maintaining the require- 
ment that the wave function be symmetric under the 
interchange of quarks 1 and 2 (recall that the other 
symmetries are handled as discussed in Sec. 



II A), the 



flavor-space part of the wave function can be constructed 
from the sum of only two components. One will be 
the product of components antisymmetric under 12 in- 
terchange in both flavor and momentum space, and the 
other the product of symmetric components. In analogy 



G 



with Eq. (4.7), the flavor-space wave function will then 



be a superposition 

(P; k x , k a ) =cos 6 D 0°*£?(P; k x , k 2 ) 

+ sin0 Z5< /, 1 vI/^(P;k 1 ,k 2 ) (4.24) 

where mi is the projection of the spin-2 spatial wave 
function in some (arbitrary) fixed direction, and the 0,1 
are the isospin operators (4.3). In parallel with (4.8), the 



spatial wave functions will be written 



mi 



k?y 2 2 m ,(k x )-krOk2) 



1 2m e » 

>k 2 r r 2 ( k ) 



sin^HFi (f) $ c (4.25) 



where ty Da is the most general L — 2 function depending 
on ki and k 2 that is antisymmetric under the interchange 
ki -H- k 2 , and ty Ds the most general symmetric function, 
and is a symmetric D-state counterpart to the S- 
state function $5. Note that the mixing between the two 
terms in ^ Da is fixed by the antisymmetry requirement, 
but that the requirement of symmetry alone cannot fix 
the relative contributions of the two terms in \E ,I,S . We 
will fix these below. 

To complete the construction of the D-state wave func- 
tion we combine the L = 2 orbital part with the total 
quark spin 3/2 wave functions. For a spin- 1/2 state with 
projection A this combination is 



Note that both of these components depend on the di- 
quark (spin-1) polarization A. 

Before moving on with our construction, we call atten- 
tion to the work of Diaz and Riska [TT], who present a 
construction of a D-state wave function along lines sim- 
ilar to ours. However, we were unable to confirm their 
results. 



2. Alternative nonrelativistic form 

To prepare for what follows, it is useful to rewrite the 
O functions in Eq. (4.29) in an alternative form. In Ap- 



pendix [Bl we show that 



0AA(ki) = 



(eX) t D U '(Mi)<re\^X), (4.30) 



where s\ is a spin-1 polarization vector of the (12) pair 
with spin projection A (already introduced above), and 
|| A) is the two-component spinor with spin projection A, 
both along the z direction, and the angular momentum 
two D matrix is 



D £m (ko = kfk™- 1 k 5 2 ^ 



3. Normalization 



(4.31) 



Dr) _ 



E 

mi—— 1 



<2m,§ M |§A) 



(4.26) 



where 77 = {a, s}, and the value of fi is fixed by the vector 
coupling coefficient, and the spin 3/2 state can be written 
as a direct product of a vector and a spin 1/2 spinor: 

|| / x>=V(lA|H|M>|lA)®||^) (4.27) 



where 1 1 A) = e\ is the spin function of the (12) pair 
(spin 1 and polarization A), and || v) is the spin of the 
spectator quark 3. It is convenient to project out the 
polarization state of the diquark, and the wave function 
we will use in applications will be constructed from the 
functions 



mi — — 1 



(2m4/x||A) k 2 l- 2m ,(k,) <1A|| M > 



4tt J2 (2m4 )U ||A>k 2 F 2rof (k i ) 

mi— — 1 

X (1 A§i>|§m) || i>). (4.28) 

The complete three-quark nonrelativistic symmetric D- 
state is then 



^(P;k x ,k 2 ) = ^(©AA(ki) - e£ A (k 2 )) $ 



D 



^(P;k l5 k 2 ) 



cos^0£ A (k) + sin^e£ A (r))$ 



(4.29) 



The wave function (4.29) depends on the momenta of 



all three of the quarks. Extracting the diquark content 
requires that we average over the internal momenta of 
the 12 pair (the diquark), just as we did in Eq. ( |4.14[ ). 

Using the new representations for the O's, it is now 
possible to separate the internal diquark variable (r) from 
the quark variable (fc). Denoting ki = k + = ^k + r and 
k 2 = k = 



±k 
2 K 



r, the D of Eq. (4.31 ) separates 



D em (k ± ) = \D em (k) + D e 



l (k,r) (4.32) 



with 



G lm (k, r) = k l r m + r"k m - U tm k ■ r (4.33) 



With this substitution, the antisymmetric component 
becomes 



^(P;k l5 k 2 



Xov, 



|A)$ D (P;ki,k 2 ). (4.34) 



Note that this component is linear in both r and k, de- 
scribing the coupling of a diquark with an internal P- 
wave, or internal angular momentum-l structure to the 
third quark in a relative P-wave, while the isovector com- 



ponent (4.29) is a sum of two D-wave components, one 



in the total diquark momentum k and the other in the 
internal diquark momentum r. 

In matrix elements, when we average over the direc- 
tion of the relative momentum, r, these components will 
give a non-zero result only when multiplied by another 
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component of the same type. Using the results from Ap- 
pendix [C] and the isospin matrix elements (4.10), the 
normalization integral becomes 



Q =eo 



1-3 



COS 2 6>D ( * 



.Da 



+e 
4e 



1-3 



sin vd 



Ds\ 



l + r 3 (3cos 2 6>£, - sin 2 M 



X'X 

S yx M D (4.35) 



with the normalization constant Md dehned in Eqs. ( C2 ) 



and (C7). Once again we observe that the correct neu- 
tron charge cannot be obtained unless Op = ir/4, imply- 
ing an equal mixture of diquark symmetric and antisym- 
metric components, just as for the S-state. The D-state 
wave function is then normalized to Hd = l/ e 0- 



4- Defining the nonrelativistic quark- diquark wave function 

The D-state presents us with some new issues in ex- 
tracting a quark-diquark wave function from the full 
three-quark wave function. Examination of Eqs. (4.29) 



and (4.34) shows that there are three orthogonal struc- 



tures, each with a different kind of average over the di- 
quark internal momentum variable r. The first of these 
accompanies the square of the 0^ A (k) term in the isovec- 
tor term. Including an extra factor of 1/4 (to simplify 
the final normalization) gives a definition 



NR 



k 4 $ 2 n = 



k.r 



d 3 k 

(2tt)3 



(4.36) 



where the integral was defined in Eq. (2.13), and to sim- 
plify the notation we let <I>d(P; ki, k2) — >• and the 
new wave function <j)r) = (/)d(P, k). Except for the factor 



of 1/4, this is precisely the same prescription used for the 
S-state wave function [see Eq. (4.14) for the arguments 
and other details], and leads to the replacement 



This function is normalized to 
d 3 k 



(27T) 



3^ 



1 

eo ' 



(4.37) 



(4.38) 



which, because of the definition (4.36), is equivalent to 
the normalization defined in Eqs. (C2) and (C7). 

The second term to be defined depends linearly on r 
and appears only in the isoscalar term (4.34). As pre- 



viously mentioned, this linear r dependence describes a 
diquark with an internal angular momentum dependent 
P-wave, which can be represented by a diquark with a 
polarization vector £„ (with v = {±, 0} the three inde- 
pendent polarization states). But this isoscalar diquark 
is described by two vectors: one due to its internal mo- 
mentum, £„, and the other due to its spin, e\. This 
term is orthogonal to all other terms, and the diquark 
content of this term can be extracted by introducing the 
following correspondence 



NR 
k.r 



k 2 rV"$ 2 5 = ( 5, m - 



NR 



k 2 r 2 $ 2 



NR 



Ar 2 

^ C P (1,-m* 
2 / ^ Wsi/ 



3 k 4 ^, 



(4.39) 



where the factor c p — l/A = 3/20 was computed in 
Appendix [C] [see Eq. (C7)] and the extra factor of 4 from 
the correspondence (4.36) has been included. Note that 



the integral over the continuous variable r is replaced by 
the sum of the complete set of polarization states £„, and 
the average over r 2 is replaced by cp k 2 . With this choice 
we obtain the correspondence 



cp 



k|C5< 



(4.40) 



with the sum over polarization states, v, to be carried 
out in the calculation of any matrix element, using the 
completeness relation 



£c£e 



Hm 7 



(4.41) 



and with the understanding that the terms linear in the 
total diquark momentum, k, will be left unchanged. The 
reason f or usi ng the momentum k (instead of r) on the 
r.h.s. of (4.40) is to eliminate all dependence on r, leav- 
ing only one function 4>d with the same normalization 
condition (4.38). 



Finally, the last term to be defined is the term depend- 
ing on e£ A (r). This term is a contraction of the spin-l 
vector, £\, of the diquark with the D-wave internal mo- 
mentum structure of the diquark. Its average can be 
represented by a new effective diquark polarization vec- 
tor, ejjAi with the familiar property (4.41). The average 
we need is 

£ / (e A ) ir D e ' m '(r)(e* A ) e D^(r)^ D 
a J k.r 



NR 



r 4 = 5 m ,, 



2c 



D 



k , r 



NR 



k.r 



d 3 k 

9"^- CDACDA J W 

A v 



8c D y—^ 1 



(4.42) 



where c 2 D = B/A = 1/16 was computed in Appendix O 



and the factor of 4 from (4.36) has again been included. 



With this choice we obtain the correspondence 



Ue* A ) e D em (r)<S> D 



V2c D 



-DA 



k 2 (/>_D 



(4.43) 



With this notation, the quark-diquark wave functions 
corresponding to the three-quark wave functions given in 
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Eqs. ( |4.29[ ) and ( |4.34[ ) become 
^i\ Q (P,fc) = C p(£X),G £m (k,C,)a m |iA) |k|^(P,fc) 

<A(^fc)=^©AA(k)0D(P,fc) 



4cd 



15 



DA a m | 2 



| A) k 2 <MP,fc), (4-44) 



where G a/3 (fc, £„) is the straightforward generalization of 
its nonrelativistic counterpart, = y— k 2 , 

®Ax(P,k) = ^iUD a P(P,k)Up(P,\), (4.49) 

and <pD(P,k) is a spherically symmetric scalar function 
of the off-shell quark momenta k^ = P — k. 



where the value of the mixing angle <f> determined in Ap- 
pendix [C] has been used. 

It is straightforward to confirm that our substitutions 
preserve the normalizations given in Appendix [C| W e 
now turn to the relativistic generalization of Eq. (4.44). 



5. 



Relativistic quark- diquark D-state wave function 



Using the ideas developed in Ref. (TUJ H2], the rela- 
tivistic analogue of the wave functions from the previous 
section arc constructed from the four-vector 



k = k 



M 2 



P. 



(4.45) 



where P is the total four momentum of the nucleon. This 
insures that k will reduce to k in the rest system. Simi- 
larly, all of the polarization four-vectors are chosen to be 
orthogonal to P, insuring that, in the nucleon rest frame, 
their time components are zero and their spatial compo- 
nents are identical to their corresponding non-relativistic 
three- vectors. In addition, we replace 



£>«' - 

G"'(k,G,)- 
D u '(k) - 
<t> D {P,k)- 

where I — > 



-g a p — - \g a p 
(4) Q 

V3U fj (P,\) 

D ap (P,k) = k a ¥ 
^ D (P,k), 



P a Pp \ 
M 2 ) 



CP 



(4.46) 



a, HI — > ft and U^{P 1 k) was defined in 
Eq. (|4.17|). With these correspondences the total D- 



state wave function of the nucleon for an on-shcll diquark 
(composed of quarks 1 and 2) and an off-shell quark (3) 
with momentum k% — P — k is, 

*aa(^ *) = ~{</>V?A k) + ^xiP k)} (4.47) 



where 



^ A Q (P, k) =V3c P (el) a G^(k, Cu)U p (P, \)\k\i/> D (P, k) 

^(PV^e^ipk^ntpk) 

- ^| e^ K U p {P, A) ~k 2 yj D (P, k), (4.48) 



6. Normalization of the relativistic wave function 

The normalization of the relativistic D-state wave func- 
tion, like the S-state, is fixed by the charge (|3.5|): 



Q = 3 E / *AAP,k)j q (Oh°*% x (P,k). (4.50) 



This condition is satisfied by the wave function (4.47), 
with components ( 4.48 ) , provided 

which is the relativistic analogue of the nonrelativistic 



1 



norm (4.38 1 



From now on we will use = \J m 2 + k 2 instead of 



E, 



D. P-state component 

The origin of the P-state component can be traced back 
to the general CST relation between the relativistic ver- 
tex function, T, and the relativistic wave function 



*(P, k) = 



1 



lq k* : 



r(p,fc) 



T(P,fc) (4.52) 



where (m q — ^ 3 ) _1 is the propagator of the off-shell quark 
with four-momentum k$ = P — k (when quarks 1 and 2 
are the on-shell spectators). In the presence of confine- 
ment, the pole at m 2 = k 2 is cancelled by a zero T, so 
the wave function can be modeled without regard to this 
singularity. For simplicity, in our previous work [5] we 
also absorbed the projection operator m q + ft 3 into our 
model of "J, but the assumption that the resulting 'J is 
a pure S-state ignores some additional structure that the 
projection operator could provide. 

Now, suppose we assume that the projection operator 
is not absorbed into the definition of the wave function. 
Then the wave function would be 

r(p,fc) 



n(p,k) 



m 2 - k 2 



(4.53) 



and the action of the quark projection operator on this 
wave function would give 



(m q + fi 3 )V$(P,k)= (m q + M- 



M — 



*UP,k) 
P-k 



M 



* s (P,fc) 



ty s (P,k) + n P %ty p (P,k), (4.54) 
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where fc was defined by Eq. (4.45), and the expressions 
were reduced using the Dirac equation, fty = sat- 
isfied by the wave function. This clearly shows how the 
Dirac projection operator leads both to a redefinition of 
the S-state and to a new P-state component. By allowing 
the mixing parameter, np, to be determined by the data, 
we free ourselves from any biases about the possible size 
of this effect. 

A P-state component could have been constructed by 
starting from either the S or D-state contributions, but in 
this paper we limit the discussion to P-states constructed 
from the dominant S-state component, as suggested by 
(4.54). This gives the following ansatz for the P-state 



*L(P, k) = % [<t>°u{P, A) - 4, 1 eTp U a (P, A)] 



x^ P (P,fc). 



The normalization of the P-state is 



(4.55) 



d 3 k 



(2^) 3 2P fc 



(-&)1> 2 P (P,k) = l, (4.56) 



where the minus sign arrises because $7°$ = — fc 2 > 0. 

Our S-state wave function is a Dirac spinor with the 
lower two components exactly zero in the nucleon rest 
frame. Conversely, the wave function (4.55) is a Dirac 



spinor with upper two components exactly zero in the 
nucleon rest frame. For this reason it has positive parity 
(the negative parity of a P-wave being cancelled by the 
negative sign from the Dirac parity operator 7 ). 

This P-wave component is of purely relativistic origin 
because it will only make nonzero contributions to matrix 
elements which have lower, relativistic components, in 
Dirac space. 



V. BREAKING ISOSPIN SYMMETRY 

In fitting the DIS cross sections, it will be necessary 
to include the possibility that the u and the d distribu- 
tions are not identical. This assumption violates isospin 
invariance. 

In the following discussion, the "u distribution" refers 
to the u distribution in protons and the d distribution 
in neutrons (we retain charge symmetry). Similarly, the 
"d distribution" refers to to the d distribution in protons 
and the u distribution in neutrons. 

To generalize the formalism to allow for this difference, 
look at the flavor wave functions (4.3 ) in more detail. Ex- 



amination of their structure shows that some of the flavor 
functions describe u distributions and others d distribu- 
tions. In particular, the isoscalar and the I = isovector 
components describe u quark distributions, since the in- 
teracting quark is always a u quark in the proton or a d 
quark in the neutron. Returning to the full three-quark 
notation of Ref . [5] , where the flavor wave functions were 
written as a direct product in the order 1,2,3, so that 



016203 — > abc, and identifying the (12) pair with the di- 
quark, so that the photon interacts with the third quark, 
it is easy to extract the separate u and d distributions 



^(ud — du) 



g (ud + du) 



proton 
neutron 

proton 
neutron 



</>V 



<t>Lox' (5.1) 



while the isovector I ^ components describe the d 
quark distribution in the proton and the u quark dis- 
tribution in the neutron: 



(uu)d proton 
(dd)u neutron 



(5.2) 



Next, we introduce different wave functions for the u 
and d quarks for each angular momentum component 



^i(P,fc)-> 



ipu (P, k) u quark 
tp^ (P, fc) d quark 



(5.3) 



where L =S, P, D. We assume that (for example) the u 
quark distribution corresponding to the spin-0 and spin- 
1 diquarks are identical, but this could be relaxed later, 
if necessary. Then, the matrix element of the charge ac- 
companying the isoscalar diquark is unchanged except for 
replacing ij)\{P,k) by [ip^(P, fc)] 2 . Ho wever , the isovec- 
tor matix element, j s -0 2 [recall Eq. (4.10)], now sepa- 



rates into two independent contributions. Allowing for 
the fact that the renormalization of the u quark charge 
(eo — > e„) and d quark charge (eo — > e° d ) might be un- 
equal, and dropping the index L, gives 

j^ 2 u (P, fc) = (r • £ )(A + §t 3 )(t • Co) el^ u {P, fc) 
= {\ + \Tz)e Jl{P,k) 



j^ 2 d (P, fc) = (r • + ±t 3 )(t • C+) e^OP, fc) 

+(r^-)a + hr 3 )(r-C)e°A 2 d (P,k) 



= (|"T 3 )e> 2 (P,fc). 



(5.4) 



Hence the generalized isovector charge operator is trans- 
formed (for any L) to 

j- fc) ^j^l(P, fc) + if tf(P, fc) 

= Ul-T 3 )\e^ 2 u (P,k)+2e° d i>l(P,k) 



+ |t 3 



e> 2 (P,fc)- e°^ 2 (P,fc) (5.5) 



If tpd = ipu (and e° = e") th en the sum of these reduces 



to the previous result ( 4.10| ), but tp^ ^ ip u the last term 
clearly gives a different result. When combined with the 
isoscalar contribution, the total is (choosing the mixing 
angle 9l = tt/A as before) 

lj A i> 2 u (P, fc) + 1 [j^l(P, k) + jf^(P k)] 
= i(2e^ + e^) + i(2 e > 2 -e> 2 )r 3 

_/ |(4e°^-«) P 



(5.6) 
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This gives the requirement that both u and d distribu- 
tions be normalized to l/e°. For arbitrary L this gives 



d 3 k 



{27r) 3 2E k 

d 3 k 
(2ir) 3 2E k 



(-k 2 ) L [^(P,k)] 2 
(-k 2 ) L [^(P,k)} 2 . 



(5.7) 



This result for the charge operator may be generalized. 
In particular, if O is any operator in isospin space, and 
the isospin one matrix elements of O are defined by 



3^0^1 = 01 

Z<t>\0<t>\ + 3 ^lO^li = 0\ 



(5.8) 



then for arbitrary angular momentum components L and 
L', we may make the replacement 



{0) L ,r^{0 + Ol)^{P,k)^{P,k) 



+0\tf\P,k)tf(P,k). 



(5.9) 



VI. SUMMARY AND OVERVIEW 



In this paper we present relativistic CST wave func- 
tion for the nucleon with S, P, and D-state components. 
These wave functions are designed to be used in calcu- 
lations where two of the quarks are non-interacting on- 
shell spectators, with the third off-shell quark interacting 
with an external probe. In such a situation the full de- 
pendence of the matrix element on the relative momen- 
tum, r, of the two on-shell spectators is contained in the 
three-quark wave function. Integrating over r will then 
lead to a new effective wave function with the two non- 
interacting quarks replaced by a quark pair (an effective 
diquark) with a mass, m s which may be treated as a 
parameter. The resulting quark-diquark wave function 
contains all of the information originally included in the 
three-quark wave function and may be used in a variety 
of calculations. 

The extraction of a quark-diquark wave function is dis- 



cussed in general terms in Sec. II B with the most general 
result for the extraction of the wave function presented in 
Eq. (2.9). In Sec. IV the various component wave func- 
tions are constructed and presented. The relativistic S- 



state quark-diquark wave function is given in Eq. (4.161 



This is identical to the S-state wave function used pre- 
viously [5J [TU] , but the extraction of this wave function 
from the full three-quark wave function, culminating in 
the nonrelativistic correspondence (4.14|, has never been 
discussed before. The relativistic P-state q uark- diquark 
wave function, also new, is given in Eq. (4.55). This 



component vanishes in the nonrelativistic limit. 

The bulk of this paper is devoted to obtaining the rel- 
ativistic D-state quark-diquark wave function given in 
Eqs. (4.47)- (4.48). This component is a superposition of 



terms with an isoscalar P-wave diquark (represented by 
the polarization vector £„), and two terms with an isovec- 
tor diquark, one with an internal S-wave structure and 



one with an internal D-wave structure. Although each of 
these diquarks has a different angular momentum struc- 
ture, they are all derived from contributions in which the 
quarks are in a relative L = 2 angular momentum state. 
Furthermore, each of these diquarks has spin-1, as re- 
quired by the coupling of overall spin 3/2 to L — 2 to 
produce a nucleon with spin- 1/2. These three compo- 
nents are orthogonal, and the only component that can 
interfere with the dominant S-state component is the lat- 
ter, which is summarized by the contribution 

*aaOP, k) -> ^ (e* A ) a D af3 (P, k)U p (P, X)^ D (P, k) 

(6.1) 



In applications when we need terms linear in the D-state, 
this is the only term that need be considered, and it can 
contribute only to matrix elements which do not allow 
for a complete integration over all directions of k (as in 
the important case of DIS). 



The normalization condition for the component (6.1) 



is 



3£/ K'y(p,k)j q (oh ^(p,k) 

A Jk 

2 



A' A 



(6.2) 



where j s was defined in Eq. (4.10) and the D-state nor- 
malization condition (4.51) has been retained. Recall- 



ing that the contribution from both the isovector D-state 
components is 1/2, we recover the previous result, also 
contained in (4.48), that the component (6.1) accounts 



for 4/5 of the total symmetric contribution. 

To allow for the possibility that the u and d quarks 
have a different angular momentum distribution, we may 
break isospin symmetry as discussed in Sec.[V] The sym- 
metry can be broken and the shapes of the u and d quark 
distributions individually adjusted as long as their nor- 
malization remains fixed. 

The wave functions derived in this paper are used in 
our discussion of DIS in the companion paper [9] , which 
provides numerical calculations indicating that the L / 
wave function components described here make impor- 
tant contributions to the proton spin. 
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Appendix A: Relativistic diquark averages 



gives 



In this Appendix we discuss the reduction of the inte- 
gral ( 2.3 ) to the form ( 2.8 1 . Begin by writing the integral 



in a manifestly covariant form and replacing the integrals 
over k\ and k 2 by integrals over the sums and differences, 
k and r, of these momenta 



(27T)' 



• [ = [ d 4 fc lC i 4 fc 2 8+{m 2 q -kl 

Jk 1 k 2 J 



)Mm? 



d 4 k J d 4 r 5+{El-(\k Q + r Q ) 2 ) 
x5 + (El-(\k -r Q f) (Al) 



Both integrals are covariant, so we will evaluate the in- 
tegral over r in the two body rest system, and use the 
two 8 functions to fix tq and |r|. This will leave the in- 
tegral over k completely unconstrained, except we define 
k 2 = s. Therefore, in the two body rest system kg = ^fs 
and the integral is 



Cirri"/ ' / </'/,- / </<->, \ - (A2) 



Ai/,'2 



As 



Finally we express the fco integration in term s of s using 



E 2 , giving the relation (2.8) 



k t k 2 



s k 



In using this formula, we must be careful to evaluate 
angular integrals involving r in the rest system of the 
two-body pair. 



Appendix B: Alternative form for the nonrelativistic 
D-wave 



To show the equivalence of (4.28) and (4.30), first in- 
troduc e the spin 3/2 projection operator [obtained by 
using (4.27) to write each spin 3/2 state in terms of a 
direct product of spin one and spin 1/2 states] 



^3 /2 = EIIm>( 



= ^(lA 2 iH|M><lAxi^|lM) 
x |lA a > <1A X |® ||i/><ii/| . 



(Bl) 



Now, starting with ( |4.30 1, insert this operator between 
the e* and D (we will show later that the contribution 
from the spin 1/2 projection operator is zero, and hence 
inserting the spin 3/2 projection operator is equivalent to 
inserting unity) , use the orthogonality of of the spin one 
polarization vectors ( 1 1 A) = e^), and insert a complete 
set of polarization vectors between the D and a. This 



liv' A' 



(B2) 

Use Eq. (2.20) from Ref. TU] (derived relativistically, but 
also true nonrelativistically, and with a sign correction) 
to evaluate the term in curly brackets 

D AuA , = (e* Al ) e D a ' '(ki)(e A 0/ 



k 2 Y 2mt (k l )(2m e lA 1 \lA'), (B3) 



(which fixes A' = me + Ai = me + fi — i>', allowing the 
replacement of the sum over A' by a sum over me) and 
use a new identity 

(|^|4,- C r|lA> = -V3(lA'i^|iA>, (B4) 

to fix v' = A — A' = A — me — fJ, + z/, which requires that 
/i = A — me and removes the sum over /i. With these 
substitutions 

X (1A| v||a*> (B5) 



(A3) where 



S miX =J2 (2m t 1 Ax 1 1 A' ) (lA'i^iA) 

v' 

x (lAiii/lf/A. (B6) 

This sum can be done using Racah coefficienttj^Jfrom 
Rose P3]) 

£<2m* lAtllA'XlA' \u'\\\)(lA l §z/|f M ) 

= R il (2m l | /x|| A). (B7) 
The Racah coefficient, from Table 1.4 on page 227, is 



Ji3 1 =>/i2W(2,l,|,|;l,f) 



5311 



4325i|2 



= 1 



(B8) 



1 Rose, Eq. (6.5a), gives 

(ji mi h mi\j') {f mi + m 2 J3 m 3 \j) (j 2 m 2 js m- A \j") 

m,2 or 7713 

= Rjiiji (ji mi j" m 2 + m-i\j) 

Here we identify ji = 2, mi = mf, j 2 = l,m 2 = Ai, j' = 
l,mi + m 2 = A', j 3 = |,m 3 = v' , j = \ , mi + m 2 + m 3 = A, 

j" = §, 7712 +W3 = /It. 
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Hence, the final answer for O is 
2 

Qfrfc) =yfa E (2m < |/x||A)k?y 2 m t (ki) 

m,f-- 1 

x (lAii/||M}||^>, (B9) 



which is precisely Eq. (4.28). 



Repeating the derivation using a spin 1/2 projec- 
tion operator in place of a spin 3/2 operator am ount s 
to replacing the state by li/i) in the sum (B7|, 

and in this case the Racah coefficient [proportional to 
W(2, 1, |, |; 1, |)] is zero, proving, as asserted above, 
that the spin 1/2 contribution is zero. 



Appendix C: Nonrelativistic averages over the 
D-state matrix elements 

Summing over the diquark polarization and averaging 
over the direction of the relative momentum of the di- 
quark, r, using the identities (where / is an arbitrary 
function) 

J d 3 rr irj /(r 2 )=^4/ rf3rr ^( r2 ) 
5>aM4)*=<5« (CI) 

A 

gives the average of the isospin-0 component of the wave 
function (14.34 1 



^E / Kv(o-MM)^(o-M^) 

\ A Jk,r 

/■INK, 

/ (\ A 'l *t E G ' m '( k ' r)G m ' m (k, r)a m | \ A) 
/ k 2 r 2 $ 2 D ( X )= *vaMj (C2) 

./fe,r 



where the nonrelativistic volume integrals were defined 
in Eq. (2.13), and we assumed t hat & p has the same 



argument x = 5k 2 + 2r 2 (see Eq. (|4. 13h) as $ s . The last 



Then, summing over the spin-1 diquark polarization vec- 
tor and using the orthogonality and normalization prop- 
erties of the O's gives 



V 1 



A', A 



/•NR. 

E / ^(05^,^)^(0;^,^) 

a Jk.r 



NR 



(cos 2 <£k 4 + sin 2 </>r 4 ) $ 2 



D ■ 



(C5) 



Now, using the ansatz (4.13) (with x — 5k 2 + 2r 2 for 



P = 0), we can find a relation between the integrals of 
k 4 , r 4 and k 2 r 2 over Integrating over r by parts first 
and then over k, gives 



POO />oo 

Nd= (k 2 r 2 ) = / k 2 dfc / r 2 dr(k 2 r 2 )$ 2 D ( X ) 
Jo Jo 



POO 

4 / 1,4 



k 4 dfc / drr b — 
Jo d X 
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5 



k 2 dfc 
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r 2 dr(r*)<l>l(x) = -{r 4 ).(C6) 



Doing the same steps, but in reverse order, gives a re- 
lation between (k 4 ) and <^k 2 r 2 ). These results can be 
summarized 



<k 4 )= A(k 2 r 2 ) = f (k 2 r 2 ) = 4A/- D 
(r 4 )=i?<k 2 r 2 ) = A( k V) = iAA D . (C7) 

Hence, both isospin averages can be written in terms of 
the single normalization constant, A/"d 



Da I 



A' : A 



5\>\N D 



A', A 



= 5y X N D (4 cos 2 



isin 2 ^). (C8) 



We choose 4> so that both normalizations are equal, which 
gives 



cos — 



line defines the D-state normalization constant. 

To evaluate the isospin- 1 average, first use the identity 

E / d%D em '(k t )D em (k t )f(^) 

= 5 m , m ~ J d 3 fc 4 k 4 /(k 2 ),(C3) 

to derive the following convenient result for the normal- 
ization of the O's 

E / <Ph&&,(ki)QZx(ki)fQ3) 

A J 

= I / d ^E <5 A 'l cT TO ^ W (k i )^ ro (k i )a m |i A) /(k 2 ) 
= <5a'A J d 3 hi k 4 /(k 2 ) . (G4) in agreement with ((C2 1. 



(C9) 



Before leaving this section, we point out that the rela- 
tions (C7) are essen tial for consistency. Using the orig- 
inal definition (4.29) for ty Da , an alternative calculation 



of the norm using the orthogonality relations ( C4 ) , gives 

- r>„ 1 2 \ 



5\>\ \ 



<k 4 > + <k 4 > 



- <5A'A \ 



k + r) 



„4 , 5i,2„2 



k 2 r 2 ) 
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